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Abstract 

The coherent conductance and current is calculated through two quantum 
dots using the Hubbard model for a single level per spin. The occurrence of 
negative differential conductance is demonstrated. The Ohmic conductance 
is calculated for dots with equally spaced levels. Transport is determined by 
matching energy levels, even when they do not occur at the charge degeneracy 
points. 
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The advance of lithographical techniques on a nanometre scale in recent years has made it 
possible to study systems that were inaccessible to experimentation before. This has opened 
the way to produce and investigate structures where the carriers are confined to one or even 
zero dimensions. This has given rise to the discovery of a number of novel effects. It has 
been shown by Reed et al. that discrete states can clearly be discerned in quantum dots, 
structures which have been confined in all three dimensions^. These technological advances 
have made it possible to study the interplay between charge quantisation effects (Coulomb 
blockade) and size quantisation effects (discrete energy levels). 

Whereas there has been a considerable amount of experimental and theoretical study on 
the conductive properties of single quantum dotsil and single and double metallic dot^, 
so far relatively little attention has been paid to the case of two quantum dots in series. 

The transport properties of a single dot are insensitive to incoherent scattering, provided 
that the broadening of the levels is smallllli. However, when two dots are connected in series, 
this no longer holds true. In this paper the coherent case will be considered where the phase- 
breaking rate is small compared to the tunnelling rate. Recently, the importance of coherence 
for transport through a single dots has been shown explicitly by direct measurement of the 
phase of the transmission coefficientli. 



I. METHOD 

The system under investigation consists of an interacting region of one or more quantum 
dots connected between leads or electron reservoirs of specified chemical potential. In general 
the reservoirs will have different chemical potentials, causing a current to flow through the 
dots. This makes it inherently a non-equilibrium process. Therefore the system will be 
best described using non-equilibrium Green's functions, now commonly referred to as the 
Keldysh method0ll3. This formalism is not only valid in the linear response regime, but also 
for higher bias voltages. A full description not only requires knowledge of the retarded and 
advanced Green's functions, as in equihbrium, but also of the 'distribution' Green's function 



2 



between the reservoir and dot sites. 



Gl^{t,t') = -iQ{t-t'){{c^{t),cl{t')}) 



(1) 



Q^(t,t') = ^e(t'-t)(K(t),4(t')}) 

G<^{t,t') = i{cl{t')c^{t)) 



(2) 



(3) 



where c' and c are the usual creation and annihilation operators. In steady state the diagonal 
elements of C — G'^ are proportional to the local density of states, whereas G*^ plays the 
role of the density matrix. They are related by the distribution matrix F, defined by 
G^ = — F(G^' — G"). In equilibrium F is simply a scalar and is identical to the Fermi-Dirac 
distribution function. 

Having introduced the non-equilibrium Green's functions, the obvious next step is to 
find an expression for the Landauer formulalll relating the current to the local properties of 
the system, such as the chemical potentials of the reservoirs, the density of states and the 
average occupation of the dots. It is assumed that the left and right electron reservoirs are 
large enough to have well-defined chemical potentials and ^r. Any interaction effects in 
the reservoirs can be neglected as a result of the screening by the free flow of charge carriers. 
Transport proceeds by electrons hopping between the reservoirs and the interacting region. 
The Landauer formula can be written a 



where {Tijjijnm = Svr X^aeL/R Pa(^)^a,n(^)V'a,m('^) matriccs coupling the interacting 
region to the reservoirs. Va,n are the hopping potentials between the reservoirs and the 
interacting region and pa is the density of states in the reservoirs. It is worth noting that 
the Green's functions in this formula are to be calculated from the full Hamiltonian including 
the reservoirs, even though the reservoirs are already present in the couplings Tl/r of the 
current equation. 

For negligible bias voltages, knowledge of the retarded (and hence the advanced) Green's 
functions allows the distribution Green's functions to be calculated. The next step is to 



2h J 



(4) 
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choose a particular Hamiltonian H and to calculate the retarded Green's functions explicitly. 
This can be done using the equation of motion method which consists of differentiating the 
Green's function with respect to time, thus creating higher order Green's functions. In 
Fourier space this amounts to the following iteration rule {h = 1): 

u;{{c;cl)) = {{c,cl}) + {{[c,H];cl)) (5) 

where c can be any combination of different creation and annihilation operators and 
^rim(^) = {{^n]cl^)) is the Fouricr transform of the retarded Green's function. The above 
equation can be applied successively to produce multiple particle Green's functions. For a 
system of two reservoirs and dots with n energy levels a multiple particle operator con- 
sists of a product of at most AinN + 2) — 1 single particle operators, i.e. the creation and 
annihilation operators for both spin orientations in each dot level or reservoir. Therefore, 
it is in principle possible to obtain a closed set of equations. Clearly, it is a daunting task 
to calculate and solve all possible simultaneous equations, so one usually resorts to making 
some approximations in order to close the set of coupled equations. 



II. SINGLE DOT 

In order to determine a realistic, yet simple, set of approximations, consider the simple case 
of a single quantum dot where only one relevant level for each spin direction will be taken 
into account. Any lower levels are assumed to be permanently occupied and higher energy 
levels are always empty. The two relevant states in the dot are labelled by the electron spin 
quantum numbers a or a. A tight binding Hamiltonian will be used with a Hubbard term 
U describing the Coulomb interaction between the spin up and spin down electrons. This is 
analogous to the Anderson model of localised impurity states in metalJi^. 

Hl= Yl (^kaclaCka + J2^aCiCa + Unana+ J2 (VkaclaCa + H .c) (6) 
kaeL,R a kaeL,R 

where and e^a are the single particle energy levels in the dot and reservoirs respectively. 
Transport between dot and reservoir proceeds by means of the hopping potential Vka- Any 
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interaction that might take place with electrons in the full lower energy levels is incorporated 
in the definition of e^. The hopping potentials V^a are taken to be real, because of the time- 
reversal symmetry. It is required that any set of approximations yields the correct solution in 
the limit of zero bandwidth, i.e. negligible hopping potential. Moreover, the single particle 
density of states integrated over all energies must equal unity. The density of states must be 
non-negative at all energies. Finally, the solution should display electron-hole symmetry. It 
should be invariant under the transformation {n 1 — n, [/ — —U, — ^ Cq, + f/, / 1 — /}. 

The crudest set of approximations (scheme 1) consistent with these condition Jll cor- 
responds to the Hartree approximation schemed where electrons of opposite spin do not 
interact directly, but only through a modified average potential field. The resulting Green's 
function is 

((r ■r^)\ = uj-e^-Ujl-jn^)) 

- e^)iLU - - f/) - So4^ - - f/(l - {n^))] ^ ^ 

with 

So.= E ^^^ = -t{TL + Tn)/2 (8) 
keL,R ^ ^fc" 

where T^r = 27rpi ^p. is the self-energy due to the tunnelling of the a electron 
and (ria) is the average occupancy of the opposite spin state. The same Green's function 
results from the method as indicated by Caroli, Combescot and Nozieres0, where the full 
Green's function is calculated as a perturbation of the system without transmission to either 
reservoir. The occupation is the integrated product of the (diagonal) distribution function 
and the local density of states in the dot. Since the Green's function is related to the density 
of states which in turn depends on the average occupation number, it means that equation |^ 
is to be solved self-consistently. It can be shown that the integrated density of states equals 
unity, subject to the condition < (n) < 1, which is obviously satisfied in any real system. 

A better approximation scheme (scheme 2) is one which neglects the simultaneous tun- 
nelling of the electron of opposite spin, but which does not decouple the higher order Green's 
functions as in approximation scheme 1. This yields the following Green's function: 



As the tunnelling of the opposite spin electron is not taken into account, its creation and 
annihilation operators are always paired together in the higher order Green's functions. 
The operators for both spins do not truly intermingle, which results in the probabilistic 
expression for the Green's function. Therefore, the same expression can be obtained by 
solving the Green's function for the isolated site and adding an imaginary part to the site 
energy to include the reservoirs, i.e. the reservoirs simply serve as a source of broadening. 

A more accurate method (scheme 3) only neglects terms involving correlation in the 
leads, valid for temperatures higher than the Kondo temperature^. The Green's function 
corresponding to these approximations isi 

{{ca, cD) 



+ N TT (10) 



The self-energies due to the tunnelling of the a electron are defined as 

Si„H= ^ + ^—Tj 1 (11) 

fcgj^^^ Vcj — + — Cfca Co" — — — (7 + eka' 

where = f{eka), — ^ ^ fi^ka) and = 1. As for scheme 2, it can be shown that 
the density of states integrates to unity for all values of (n). 

Since the couplings to the reservoirs differ by a constant factor the distribution matrix 
reduces to a diagonal matrix. 

^ ^ fLTL + fRTn ( ^ 0\ .^2s 
Tl + Th [01) ^ ^ 

As there is only one level for each spin the current expression ^ takes an especially simple 
form where the current is proportional to the density of states between the Fermi levels of 
the reservoirs. 



The Ohmic conductance is given by the integrated product of the density of states and the 
derivative of the Fermi function. At low temperatures the derivative of the Fermi function 
approximates a 5-function. A Lorentzian density of states with its broadening caused solely 
by tunnelling to and from the reservoirs thus gives rise to a conductance peak of height e^/ h, 
provided that the dot is symmetrically coupled to the reservoirs. 

The three approximation schemes are compared in figures |l]- ^ for the case where the 
temperature equals the broadening due to the reservoirs. In order to distinguish between 
the electrons in the dot, a magnetic field is applied giving rise to a Zeeman splitting of the 
energy levels (ea = 0.0?7, = O.lf/). 

Figure || shows the current and the associated differential conductance through the dot 
when a bias voltage is applied across it, keeping the chemical potential in the right reservoir 
fixed {hr = —0.5U). All three approximation schemes yield qualitatively the same results, 
although scheme 3 (and scheme 2 to a lesser extent) has a larger associated broadening. 
This is a direct consequence of the fact that scheme 3 takes the tunnelling of the oppositely 
charged electron into account. 

In figure ^ the Ohmic conductance is plotted. The approximation schemes give identical 
solutions for the shape and height of the major conductance peaks, but especially approxi- 
mation scheme 1 gives a different result for the amplitude of the minor peaks. Moreover, an 
investigation of the occupation plot |^ shows that for scheme 1 the average occupation of the 
higher energy spin state can exceed that of the other spin state for intermediate values of 
the chemical potential of the reservoirs. This is clearly unphysical. Both these observations 
justify the conclusion that scheme 1 should be dismissed as a good set of approximations. 

As expected the conductance plot shows peaks at the single particle energy levels and at 
the same energies displaced by U. It is clear that two peaks are suppressed, which can be 
explained as follows. The first peak occurs when transport takes place through the lowest 
level ea- As the chemical potential in the reservoirs is raised, the state will start to fill up. 
When the chemical potential lines up with e^, level will be mostly occupied (see figure 
^ thus putting the Coulomb blockade in place and impeding the flow of electrons through 
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level ea- This results in a highly suppressed conductance peak. Because of the electron-hole 
symmetry a similar situation occurs at much higher values of the chemical potential. 

Approximation schemes 2 and 3 seem to give results which are closely in accordance with 
each other. The main difference is that approximation scheme 3 produces a larger associated 
broadening. Although slightly better, scheme 3 is harder to manipulate numerically, whereas 
scheme 2 is much more stable. Therefore approximation scheme 2 will be employed here. 

III. TWO DOTS 

The single dot Hamiltonian can easily be extended to describe two quantum dots in series. 
As before, assume that there is only one spin-split level per dot and that the electrons 
experience an on-site interaction energy C/i/2. 

ka£L,R « 

+ Uiniania + U2n2an2a + Yl ^M{c\o,C2a + H.C.) 

a 

+ E (^fea4aCla + H.C.) + (^fea4aC2a + H.C.) (14) 

kaEL ka&R 

As before the dots are coupled to the reservoirs by the hopping potentials Vka-i but there is 
also some coupling Vm between the dots. In this representation the coupling matrices for 
each spin are given by 

Tl = ; r« = (15) 

V 0/ Vo Vr) 

Now that there is more than one dot in the interacting region, it is clear that the correlations 
between the occupation numbers of the dots become important, since the occupation of one 
dot will depend on the occupation of the other. The equivalent of approximation scheme 1 
actually ignores this correlation. Here the approximation scheme 2 will be employed, where 
the reservoirs are included in the form of a self-energy. The single particle Green's function 
matrix for spin a is given by 
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where 




— |r/j. As for the single dot case, the Green's functions have 



a straightforward probabihstic interpretation. They treat the electrons of opposite spin a as 
static entities whose main influence is to change the effective site energy of the electrons of 
spin a by V . Each Green's function is simply a sum over effective non-interacting Green's 
functions each weighted by the probability of a particular realisation of the occupation 
numbers for the opposite spin states. Of course the single particle energy levels have to be 
adjusted to correspond to the correct occupations of the opposite spin states. 

This form of the Green's function has some important consequences. 
• First of all the distribution, retarded and advanced self-energies can be written in terms 
of the coupling matrices: = i^jjXh + Jb^r) and T7 — S° = —i{TL + T/j). Substitution 
of the equalities G< = G'^S<G'' and - G° = G''(S'^ - S'^)G° into equation | yields a 
condensed form of the current formula for non-interacting systems. 



• Secondly, the coupling matrices are no longer proportional, which results in a non-diagonal 
form of the distribution matrix at finite bias. 




(17) 



Using equations for the double dot this reduces to 




a 



(18) 




(19) 
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• Thirdly, the probabihstic form of the Green's functions warrants a statistical approach for 
the multiplication of Green's functions. Since the Green's functions are probabilistic sums 
over non-interacting Green's functions, products of Green's functions must be expressed as 
probabilistic sums over the products of the appropriate non-interacting Green's functions. 

• Fourthly, higher order Green's functions, needed to calculate correlations of the type 
{nian2a)i can also be reduced to a sum over non-interacting Green's functions. For a non- 
interacting system it can be shown that 



In the presence of interaction on the dots, the correlation {nian2a) is obtained by taking the 
probabilistic sum over all possible occupation realisations of the opposite spin states. 

Taking the four points above into account, the current and conductance characteristics 
can be calculated. Figure ^ shows the Ohmic conductance through two dots which are 
identical in all respects apart from a relative energy off-set of OAU. In order to be able to 
interpret the conductance plots more easily the average occupation numbers for the four 
single-particle states are also plotted. 

In the limit of negligible coupling between the dots, the occupation for each dot is 
determined solely by the coupling to the adjoining reservoir. This is identical to the single 
dot case ( c/ figure H). The Ohmic conductance displays peaks at all the single particle energy 
levels and at the same energies displaced by U. Similar to the conductance through a single 
dot (see figure H), the peaks come in pairs separated by the Zeeman energy — e^. As 
it has been explained for the single dot, one of the peaks of each pair is suppressed as a 
result of the Coulomb blockade in the dot. However, it is apparent from figure ^ that there 
is also a modulation of the height of pairs of peaks. For instance the first major peak is 
noticeably larger than the second peak. The current path that gives rise to the first major 
peak encompasses the levels eia and e2a- At the chemical potential at which the second 
major peak occurs, the first dot will be virtually completely occupied and the Coulomb 
blockade will be in place. Therefore the dominant contribution to the current will come 




(20) 
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from the energy levels eia + U and e2a- In the case of figure ^ the dominant levels for the 
second peak are energetically further apart than those for the first peak, which explains why 
the second pair of peaks is slightly suppressed. In short, the form of the conductance plot 
for negligible inter-dot coupling Vm can be understood by realising that electrons contribute 
to the transport only when the Coulomb barrier is overcome in both the first and the second 
dot. 

As the coupling between the dots is increased, the energy levels in the dots start to 
interact, repelling one another. The conductance peaks are shifted in energy. When the 
coupling between the dots becomes appreciable compared to the energy level difference 
^2a — ^la for the decoupled dots, the number of peaks in the conductance can increase to a 
number that is greater than the total number of single particle energy levels. This is a result 
of the fact that the amount by which a level is repelled depends on the energy difference 
with the level by which it is repelled. For instance, an original level at an energy eia can 
be repelled to two different resultant energies depending on whether it interacts with level 

or + U . This causes a conductance peak both to shift in energy and to split. This is 
particularly significant when a level is repelled into opposite directions by the two possible 
levels it could interact with. This happens at intermediate values of the chemical potential. 
This explains why splitting can start to be observed for the middle peak pairs of figure ^ 
but not for the peaks at more extreme values of the chemical potential. 

In figure |^ the current through the two dots is plotted as a function of the chemical 
potential in the left reservoir, keeping the chemical potential in the right reservoir fixed. 
The I-V characteristics are studied as the coupling between the dots is increased. 

In the limit of negligible coupling the occupation of each dot is completely determined 
by the tunnelling of electrons to and from the adjoining reservoir. The barrier between the 
dots is clearly the current limiting segment, so that the current is proportional to V^,^ (see 
figure |a). Transport can proceed because the wavefunctions corresponding to the levels 
in the first dot can leak slightly into the second dot and vice versa, thus creating current 
paths. This explains why there is a marked increase in the current at /il = eia, /^l = ^2a and 
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IJ'L = ^2a- Note that no current step occurs at /i/, = ei^, even though an extra current path 
becomes available. This is most easily comprehended by drawing a distinction between the 
number of current paths and the number of effective current paths. A current path consists 
of a composite state which is a mixture of a state in the first dot and one in the second dot. 
A current path can be said to be effective when it can be used constantly. For example, it 
is true that a second current path becomes available at /il = ei^. However, the number of 
effective paths remains fixed at 1, since the Coulomb blockade in the first dot prevents the 
paths from being used simultaneously. 

The most notable feature of figure §a is the region of negative differential conductance 
around /i/, = eiQ, + V . This is the energy at which the single particle levels in the first dot are 
both likely to be occupied. The dominant current contribution arises from the interaction 
of the levels ei + t/ with the levels €2. As these states are energetically further separated 
than the levels ei and 62, this means that there is less overlap of the wavefunctions in the 
two dots which accounts for the drop in current. 

Finally, the more realistic case of the Ohmic conductance through a double dot with an 
infinite number of energy levels will be calculated. The simple case of negligible broadening 
and small inter-dot coupling will be considered. The coupling to the reservoirs in the 'site' 
representation is given by the following matrices: 



























(21) 



Rj 



The matrix is subdivided into four infinite sub-matrices, corresponding to coupling either 
between levels of the same dot or between levels of different dots, e.g. the top left sub- 
matrix includes all couplings between levels of the first dot. The sub-matrices Fl//? are 
defined as a matrices whose elements are all given by Fj^//j. As before the current will be 



determined using equation 0. Substitution of the coupling matrices and taking the limit of 
an infinitesimal (Ohmic) bias across the system yields the conductance equation. 
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r-^ f.1 ^^^^ 
^~ hJ """^AkBT cosh' 



^2kBTl 



iGdotl j6dot2 



(22) 



Analogous to expression |TB| the conductance only depends on the Green's functions between 
states in different dots. Treating the inter-dot coupling Vm as a small perturbation allows 
these off-diagonal Green's functions to be written as Gl^iu) = G[ (07)^4/(7^(^7) so that 



G 



h 



dcu 



r2 

'M 



E Gliuj] 
iedoti 



E G^c^) 



jedot2 



(23) 



In this approximation the diagonal elements Gl{uj) of the Green's function are identical to 
those of the single dot with multiple level^. Hence 
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E Gliuj) 




igdotl 





E E' 



(24) 



where Pi{Ni) is the probability that the first dot contains Ni electrons on levels other than 



ej. In the limit of negligible broadening this becomes 

2 



E Gri^) 

iedoti 



2^ E E^^(^i)^(^-e*-^if^i^ 

jgdotl Ni 



+ r. E E 



viedotl A^i 



(25) 



The first part of the expression shows the behaviour at energies coinciding with an energy 
level whereas the second part gives the off-resonance contribution. A similar expression 
may be obtained for the Green's functions of the second dot. Note that the limit of small 
broadening excludes the possibility of energy levels in the two dots matching up exactly 
(this case will be examined later in more detail), so that the conductance formula ^can be 
rewritten as 

2Wl, 



X 



iedoti Ni 



+ Y: T.^LPAN2)S{uj-e,-N2U2) 

jedot2 N2 



E E 



LlJ - e,; 



(26) 
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In figure ^ the Ohmic conductance has been plotted for a system where the dots have a 
different level spacing and Coulomb interaction. As expected at low temperatures, the con- 
ductance peaks are largest at values of the chemical potential where the charge degeneracy 
points of the two dots are close in energy. At the charge degeneracy point between the oc- 
cupations N and + 1 the (A^ + l)th single particle level is dominant. As the temperature 
is raised there is a higher probability of other levels contributing to the current. This has 
a particularly strong effect when any of the newly available subsidiary levels very nearly 
match up (as at ;U ~ —A.7U). Otherwise a rise in temperature simply causes the peaks to 
be smeared. At higher temperatures the periodicity of the smaller dot U2 + A^j can easily 
be estimated from the spacing of the peaks in the conductance plot. 

At this stage it is useful to analyse the height and shape of the conductance peaks in more 
detail. Moreover, in order to compare the coherent case with other tunnelling mechanisms, it 
would be helpful to obtain an expression for the current through the double dot at finite bias 
as a function of the gate potential of one of the dots. At low temperatures these quantities 
can be investigated by just considering the dominant levels in each of the dots. Their energy 
difference is simply the separation A between the nearest charge degeneracy points. Since 



can 



this is now a non-interacting system, the Green's functions G12 of the current formula ITB 
be calculated exactly. It is also possible to include the effect of non-negligible dot-reservoir 
and inter-dot coupling. This yields 



J = — du! 



1 2 



(27) 



hJ [^2 _ + AVi, + TLTn)/^f + i [uj{Tl + Tn) + ^Tl - Tr) 
First, consider the peak characteristics of the Ohmic conductance. This is obtained simply 
by differentiating with respect to fii and setting fJ^L = fJ^R = Thus, for = Tr = T, one 
obtains 

hi 4fcsTcosh2(i^) [4a;2- A2-4V,^ + r2]2 + 4r2(A2 + 4F|f) ^ 
This shows that the conductance is determined by the product of a thermal broadening factor 



and a term which is a Lorentzian in uP' of width FA where A = \J /S? + 4V^j. Depending on 
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the relative width of the two terms, two hmiting cases will be of relevance. 

In the limiting case of small reservoir-dot coupling /c^T y/AT the conductance can be 
rewritten as 

^ fcosh-^(ll^) ^ ifr>A 

2ksnA^ + 4v- + P) 1 r^^^vf?^) if r < A ^ ^ 



2 I 2fcBT 

which is clearly consistent with expression ^ for negligible reservoir-dot and inter-dot cou- 
pling. This shows that, as a result of the level repulsion and the broadening, the conductance 
does in fact not diverge at A = in this limit. Moreover, the peak height is inversely pro- 
portional to the temperature, as was the case for the Ohmic conductance through a single 
dot. 

The question whether the conductance diverges at all will now be considered by taking 
the opposite limit of vanishing temperature ksT <^ VAF. It appears that the conductance 
either has a single peak or a split peak. 



at/i = ifr>A 
at/i = ±iJA2 + 4V;2 -P ifr<A 



,30, 



A2 + Wi, 

This shows that the conductance will not diverge but has a maximum value imposed by the 
conductance quantum e^//i, which can only be reached when A = 0,rL = r/i; = r and Vm > 
r/2. For asymmetric tunnelling barriers to the reservoirs, the conductance is maximised at 



Vm = \\J^'l + and gives a maximum conductance of {e^ / h)'Sr lT pi{V\ + T']^/{Ti + Vji)^. 

Secondly, an expression will be derived describing the current peaks which occur at finite 
bias when the gate voltage of one of the dots is raised. As before, only a single level per dot 
will be taken into account. Assuming that the energy window — fJ'R is sufficiently large. 



the integral in the current equation 27 can be solved analytically by finding the residues of 



the poles of the integrand in the Argand plane. This yields 

w A2 + {Tl + r^) V4 + vljiTL + Th) vr^iH ^ ^ 

The current peak has a Lorentzian line shape with a width that is at least as large as the 
combined width of the single particle levels in the two dots. This is consistent with the 
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experimentally observed line shapeta. In the limit of large coupling the same maximum 
current results as for the single dot case. 

IV. SUMMARY 

In conclusion, a coherent system of one or more weakly coupled quantum dots with Coulomb 
interaction can be well described using Green's functions of a probabilistic nature. In a 
double dot, the modes of the individual dots are coupled together. The resulting level 
repulsion can lead to a complex form of the Ohmic conductance and the current. Regions of 
negative differential conductance are likely to occur. The Ohmic conductance through two 
dots with multiple levels is dominated by the matching energy levels. 

The maximum conductance through a double dot /h is reached at low temperatures in 
the symmetric case of aligned levels A = 0, identical coupling to the reservoirs V l = V r = V 
and Vm > T/2. 

The peaks that occur in the current as the gate potential of one of the dots is varied 
while keeping the voltage bias across the dots fixed have a Lorentzian line shape with a 
width at least as large as the combined broadening of the individual dot levels. 
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FIG. 1. Current and associated differential conductance through a dot with a spin-spht level 
for the various approximation schemes (/Xij = —0.5U, ksT = Vl = Tr = O.Olt/, = 0.0?7, 
= O.IC/, 7o = 2e^TLTR/n{TL + r,^)) 
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FIG. 2. Ohmic conductance through a dot with a spin-split level for the various approximation 
schemes {ksT = Vl = Tr = O.OIC/, = O.OU, = O.IU). 



18 





A 


— scheme 1 






/ t 


scheme 2 




- 




scheme 3 






i' ("^ 
















\ 




y 1 1 


II 

: ij 
1 


/ 

// 
If' 




V ^ 



-0.5 0.0 0.5 1.0 1.5 

chemical potential (U) 
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{ksT = Fl = Tr = O.OlU, Ea = 0.0J7, = Q.IU). The two sets of curves refer to the occupation 
of the a and a states. 



19 



(a) 



Vm = 























■J 




J'' , , , , V' 







1.0 



0.5 




0.0 0.5 1.0 1.5 2.0 -0.5 0.0 0.5 1.0 1.5 
chemical potential (U) chemical potential (U) 

FIG. 4. Ohmic conductance through two dots each with a level for both spins, using various 
coupling constants Vm {Ui = U2 = U, ksT = Tl = Tr = O.OIC/, eia = 0.0J7, eia = O.IC/, 
^2a = 0.4[/, e2a = 0.5C/). In the occupation plots the solid lines correspond to the states of spin a, 
the dotted lines to the states of spin a. 
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FIG. 5. Current through two dots each with a level for both spins, using various coupling 
constants Vm {i^r = -2.0U, Ui = U2 = U, ksT = Tl = Tr = O.OIC/, eia = O.OU, da = 0.1[/, 
^2a = 0.4[/, e2a = 0.5C/). In the occupation plots the solid lines correspond to the states of spin a, 
the dotted lines to the states of spin a. 



21 



6.0 



o ' o ' •o' d d» 6 6 •b 'o '• ' o ' o ' o 



^ d« 6 6 « o 



CD 
O 

« 

I 2.0 

c 
o 
o 



0.0 



kT = 0.01 U 
kT = 0.02 U 
kT = 0.05 U 
kT = 0.10U 



\ 

\ 



-4 -2 2 4 6 

chemical potential (U) 

FIG. 6. Ohmic conductance through two dots with multiple levels in the limit of small coupling 
r,VM {Ui = OAIU, U2 = U, Aei = 0.13C/, = 0.32U, = Tr = F). The empty and filled 
circles indicate the positions at which the average occupation increases by one for dot 1 and 2 
respectively. 
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